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Abstract 

> 

^vq , The role of Aharonov-Bohm effect in quantum tunneling is examined when a po- 

tential is defined on the S 1 and has A-fold symmetry. We show that the low-lying 
energy levels split from the TV-fold degenerate ground state oscillate as a function 
of the Aharonov-Bohm phase, from which general degeneracy conditions depending 
on the magnetic flux is obtained. We apply these results to the spin tunneling in a 
spin system with iV-fold rotational symmetry around a hard axis. 
03.65.Bz, 03.65.Sq, 75.45.+j, 75.10Dg. 



1 Introduction 



Quantum tunneling in mesoscopic spin systems with large spin has attracted 
much attention over a decade, because of its fundamental and practical inter- 
ests[l]. From a fundamental point of view, these systems reveal many inter- 
esting quantum mechanical phenomena such as the tunneling of a large spin 
out of a metastable potential minimum, termed as the macroscopic quan- 
tum tunneling, and coherent spin tunneling between classically degenerate 
potential minima, the macroscopic quantum coherence[2]. A more surprising 
phenomenon in these systems is the perfect quenching of the tunnel splitting 
between two degenerate ground states. This phenomenon was predicted the- 
oretically [3, 4] and observed in recent experiment [5] in the spin system which 
has two-fold degenerate classical ground states. 
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The quenching of spin tunneling in a spin system is based on the quantum 
interference of the topological phase, known as Berry phase. More explicitly, 
consider a biaxial spin system which has two-fold symmetry about a hard axis 
so that there are two degenerate ground states along each direction of an easy 
axis. In the semiclassical approach, the topological phase formed by two sym- 
metric instanton paths which connect the two degenerate minima in opposite 
directions gives rise to a destructive interferenece, and hence leads to a van- 
ishing tunnel splitting for half-integer spin (i.e., the Kramers degeneracy) [3] . 
When a magnetic field is applied along the hard axis the two instanton paths 
are still symmetric and wind the hard axis in opposite directions to give the 
topological phase interference which depends on the applied field. The tunnel 
splitting in this case oscillates as the field varies, which is responsible for the 
field dependent quenching effect [4]. 

In this letter we show that the quenching of the spin tunneling can be under- 
stood from the viewpoint of Aharonov-Bohm(AB) effect [6] when the effective 
potential for the spin tunneling can be considered as a two-fold symmetric 
potential on S 1 . More generally, when a potential is defined on the S* 1 and 
has an A-fold symmetry, we investigate the effect of the AB-phase on level 
splittings from the A-fold degenerate ground state due to the tunneling. Using 
the semiclassical methods (i.e., the instanton approach), we develop an AB 
tunneling amplitude for the tunneling of a particle in the iV-fold symmetric 
potential on the S 1 , from which the level splittings are found to oscillate as 
a function of the AB-phase. Furthermore, general degeneracy conditions de- 
pending on the AB-phase are also found from the oscillating level splittings. 
By applying these results to the spin tunneling in a spin system with A-fold 
rotational symmetry around the hard anisotropy axis we derive general degen- 
eracy conditions depending on both the magnetic field and spin. When N = 2 
our results reproduce the same quenching conditions found in references [3,4], 
which confirms that the Aharonov-Bohm effect leads to the quenching of the 
spin tunneling. In the following Section, starting from a one-dimensional pe- 
riodic potential, we construct the AB tunneling amplitude in S 1 . The energy 
eigenvalues depending on the AB-phase and degeneracy conditions will be de- 
rived from this amplitude. In Section 3, we present the application of these 
results to the spin tunneling problem. Finally, there will be a summary in 
Section 4. 



2 AB tunneling amplitude in S* 1 

Consider a particle in a one-dimensional periodic potential V(0) with period 
2n and global minima at <f) = 0, ±2tt, ±47r, • • •, so that the Hamiltonian is 
7i = + V(<f>), where p = —id/dcj) is the momentum operator, and M is the 
particle mass. The transition amplitude for a tunneling between two adjacent 
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potential minima can be obtained from Euclidean time formalism. Introducing 
an Euclidean time r = —it, the Euclidean amplitude for the tunneling from 

(j) = at t' = —r/2 to = 2n at r' = r/2 is expressed as 

/C pp (2tt, r) = (2vr | e-^ | 0> pp = J V<j>(r') exp [-^[0(r')]] , (1) 



where "PP" stands for the periodic potential, and 

r/2 



-r/2 L V 7 



(2) 



is the Euclidean action for the periodic potential V(0). Since we are interested 
in the tunneling between two neighboring ground states the limit r — > oo is 
assumed. In the semiclassical approximation lC pp {2n, r) can be calculated 
using the dilute gas approximation and expressed by 

/C PP (27T, r) = J2 <*ni-n 2 ,lC?ni,n 2 (27r,T), (3) 



where ^ ni n2 (27r, r) is the contribution of n\ instantons and n 2 anti-instantons 
to the Euclidean amplitude. The general form of this is given by [7] 



p -(n 1 +n 2 )S c i 1—7 

'ni,n2 , 

n!!n 2 ! V 7r 



(27r,r) = - ,„ , (Drr^e-^J^ (4) 



where is the Euclidean classical instanton action, u is equal to 

evaluated at the potential minima = or 2n, and D is some real number 

which arises from the collective coordinate treatment of time translational 

symmetry [8]. The explicit form of D, however, is not necessary for further 

consideration. 

Now, using the identity 

S m -n 2 ,N = / g e -*(n™-^ (5 ) 



with N = 1, we can show that 



- ( I - 2Z>e~ Sc, cos0 ] r 



(6) 
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Since runs from — oo to oo the lowest energy levels are infinitely degenerate 
when the tunneling is absent. However, when the tunneling exists, these levels 
will be split and form a continuous band of energy levels. Introducing 7i \ 
ipd >— E$ | ^ >, where E# is an energy level in the lowest energy band, and 
| > is the corresponding eigenfunction, the Euclidean transition amplitude 
for the tunneling, in the limit r — * oo, can also be written as 

IT 

/C pp (2tt, t) = J d& < 2n \ >< \ > e~ E » T . (7) 



Comparing this with Eq. (6) we can see that the spectrum of the lowest energy 
band is 

E# = - - 2De- s *cos$, (8) 

and the wavefunction ip^((f)) satisfies the relation 

^(2tt) = e*%(0). (9) 

As we can anticipate from the periodic property of the potential V(4>), Eqs. (8) 
and (9) represent the well-known band structure and Bloch theorem for an 
electron in one-dimensional periodic lattice. 

If we impose a periodic boundary condition on the wavefunction such that 

^(2tt) = MO) (10) 



the above potential can be considered as a particle moving on a circle. In 
this case, since is restricted in the range < < 2ti, the periodic potential 
V{4>) is reduced to a 1-fold symmetric potential on the S l geometry. When the 
periodic boundary condition (10) holds the phase factor $ in Eq. (9) becomes 
zero, so that we can insert into the integral in Eq. (6). The Euclidean 
amplitude for the tunneling in the S 1 geometry is then obtained to be 



IC S (2tt,t) 



1 

2n V 7r 



cxp 



- 2De~ Scl ) r 



11^ 



where we have changed the superscript from "PP" to "S* 1 " to denote the 
S* 1 geometry. Notice that, since there are no degenerate levels in the 1-fold 
symmetric potential, the Euclidean amplitude provides only the lowest energy 
level E — | — 2De~ Scl instead of giving an energy band. 
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We now consider the effect of the Aharonov-Bohm phase on the tunneling of 
a particle in the 1-fold symmetric potential on the S 1 geometry. To this end 
we assume that an external magnetic field H is applied along the axis of the 
S 1 . The wavefunction then acquires a phase factor $ (i.e., the magnetic flux) 
per winding due to the AB effect. In this case, the ordinary winding number 
representation of the Euclidean amplitude can be constructed by inserting 
Z)m=-oo — m<&) to the integrand in Eq. (6). In the problem of tunneling, 
however, we need only one phase factor $, which means only m — 1 sec- 
tor contributes to the transition amplitude for the tunneling. The Euclidean 
amplitude for this case, which we call AB tunneling amplitude, is then 

(| - 2De~ SA cos $J t , (12) 



/C AB (27r,r) = — ^/-exp 



where we have inserted a subscript "AB" to specify the inclusion of the 
Aharonove-Bohm effect. From this it can be seen that the energy eigenvalue 
becomes 

E($) = ^ -2De- Sd cos$ (13) 



which is a periodic function of the flux $ with period $ = 2tt. This implies that 
the energy level oscillates by varying the magnetic field (i.e., the AB-phase). 
For a given value of $ the energy level is shifted up by 2De~ SrA (l — cos $), and 
the shift becomes zero whenever <3> = 2nn, where n is a non-negative integer. 

Let us generalize the above argument to the case of iV-fold (N > 2) symmet- 
ric potential V/v(0) i n the S 1 geometry. The corresponding one-dimensional 
periodic potential has a period of 2n/N and minima at = 0, ±2-ir/N, ■ ■ ■. In 
this case, the transition amplitude for the tunneling between adjacent minima 
is represented by JC pp (2ti/N,t), and the contribution of rt\ instantons and 
U2 anti-instantons is Q ni ^ n2 (2n/N,T) which has the same form as in Eq. (4). 
Our task is to find a general AB tunneling amplitude /Cf B (2ir/N, r). Before 
we do this, let us first calculate /Cf B (27r, r) to see how the instanton and anti- 
instanton contribute to the AB tunneling amplitude for the iV-fold symmetric 
potential. Following the same analysis as in the 1-fold case we can calculate 
this amplitude within the dilute gas approximation. The only difference from 
the previous case is to replace 5 ni _ n2; i in Eq. (3) by 5 ni - n2 ,iv- Using the iden- 
tity in Eq. (5) and inserting 8{d — <&/N) into the integrand in Eq. (6) we can 
obtain 



^ n l5 n 2 N 

r , , > s- ( i' 

exp — ( — — 2De brA cos 



2tt 



N 



(14) 
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This shows that, through the tunneling from = to = 2tt in the TV-fold 
symmetric potential, the instanton and anti-instanton behave as if they carry 
phase factors ^ and — respectively, up to overall phase factor in spite of 
their classical nature. 

We are now in position to find the AB tunneling amplitude for the TV-fold 
symmetric potential. The observation mentioned above makes it possible to 
compute this quantity without overall phase factor which is irrelevant for the 
present problem, especially for the computation of the level splitting due to 
the tunneling. Since the instanton and anti-instanton carry their own phases, 
the multi-instanton contribution to the AB tunneling amplitude for the iV-fold 
symmetric potential can be described by 

k ab (|, r) = E £ PZ„e K „ (|, r) e-*<»™> (15) 

V i V / ni T12 V i V / 



where Pn im represents all possible sequence of instantons and anti-instantons. 
For N >2 this can be expressed as 



= E 



n,m=0 



N-1 



S ni ,Nm^n2,Nn+N-l + ^n 1 ,Nm+l^n2,Nn + E ^n 1 ,Nm+'q5n2,Nn+'q-l ,(16) 

r?>2 



where the summation over 77 exists only when N > 3. Using this relation and 
the expression of G m>n2 (2ir/N, r) in Eq. (4) we can write the Eq. (15) as 



AT-1 



fN,o(xi)fN,N-l(x 2 ) + fN,l(x 1 )f N fi(x2) + E /^(^O/jV^-l (^ 2 ) , (17) 

r?>2 



where 



Xi = £9 = -De 51:1 f cos — + 
2 V TV 



' sm — ) 7 . 



(18) 



and 



OO x Nl+cr 



(19) 



which satisfies following relations; 

d N -°f N ,o(x) 



InA x ) 



dx N -° 



d N f Nfi (x) 
dx N 
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The boundary conditions for the differential equations are given by 



/iv,o(0) = 1, f N ,*>i(0) = 0. 



(21) 



Using these boundary conditions we can solve the differential equation and 
find 



_2_ 

N 



N 



cosh 2; + cos sm 77-) ex P (% cos ^r) 

N-l 

e x + 2 c °s sin ^fr) ex P ( x cos ^r) 



, if iV=even; 



if N=odd. 



f22) 



Substituting these solutions into the Eq. (17) and performing the summation 
over rj for each case of iV we obtain 



i (2tx \ 1 fio sr^ 
AB ( — > r ) = — l/lL eX P 



N 



NV 7T 



it 



^ o $-27r/c\ .27rA; 

2De~ bcl cos r - % 

2 N N 



.(23) 



Thus, the energy eigenvalues due to the tunneling in the TV-fold symmetric 
potential on the S^-geometry are 



(24) 



where fc is an integer allowed in the range 



iV , 

"77 < k < TT- 
2 ~ 2 



(25) 



Unlike the 1-fold symmetric potential, there are now A" eigenvalues for a given 
N, of which all depend on the AB phase $. The reason for the existence of the 
several eigenvalues is obvious from the fact that the potential Vn(4>) has A" 
minima within 2n, and hence there are Af-fold degenerate lowest levels which 
are split by the tunneling. Note, however, that, since each level oscillates as 
a function of the AB-phase, the levels may intersect each other to become 
degenerate states at certain values of $. To see how this happens we consider 
the difference between any two levels corresponding to k and k' {k ^ k'), 
respectively: 



AE 



N : 



„ o $-(k + k')ir 

■■ 4De~ Scl sin K — '— 

N 



sin 



ir( k - k') 
~N ' 



(26) 
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From this and the Eq. (25) we observe following general features for the de- 
generacy depending on $: 

Y double degeneracies at $ = (2n + l)7r 

> when N = even; 

Y — 1 double degeneracies at $ = 2nn ] (27) 

double degeneracies at $ = mr, when A = odd, 

where n is a non- negative integer. Note that, in all cases, the degree of the 
degeneracy is at most double. This implies that a complete suppression of the 
tunnel splitting (i.e., the quenching) is possible only when N = 2. From above 
conditions, this quenching effect occurs whenever the AB-phase becomes an 
odd-integer multiple of it. 



3 Application to spin tunneling 

As an application of above results, let us now consider the spin tunneling in a 
spin system which has an A-fold symmetry about a hard axis. Especially, we 
will look into the cases of N = 2,3,4,6 which can be realized in orthorhom- 
bic, trigonal, tetragonal, and hexagonal structures, respectively. We choose z 
as the hard axis so that xy is the easy plane. We also assume that an exter- 
nal magnetic field is applied along the hard axis. With this choice the spin 
Hamiltonian which displays the A-fold symmetry can be written as 

H N = AS 2 Z ~ C(S^ + S N )- gfi B H z S z , (28) 

where N — 2, 3, 4, 6, g — 2, /i# is the Bohr magneton, and A > C > 0. For 
simplicity, we have neglected the higher order uniaxial anisotropy terms S z 
and . To study the spin tunneling within the semiclassical methods we use 
the spin coherent state path integral approach. 

In the spin coherent state representation the anisotropy energy corresponding 
to the Hamiltonian (28) is given by 

g N (6,(j)) = < n \n N \n> 

= AS 2 (cos 2 6 - \S N - 2 sin^ 6 cos N<f) - h cos 0) , (29) 

where | Q, >—\ 8,<f) > is the spin coherent state with 8, <p being the polar and 
azimuthal angles, respectively, A = 2C/A, h = gfisH z /AS, and S is the spin 
which is assumed to be large so that the semiclassical methods can be applied. 
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Since we have chosen the z axis as the hard axis XS N 2 < 1 should also be 
assumed. The energy £n(6, <fi) exhibits TV-fold degenerate classical minima at 
= e Q , <p = 0, 2n/N,~; 2(N-l)n/N, where 9 = ir/2 for h = 0, and decreases 
smoothly to as h is increased. Since is defined in the range < < 2ir, the 
spin tunneling from (9 ,0) to (9 ,2n/N) is essentially same as the tunneling 
in the iV-fold symmetric potential on the S 1 geometry. 

To see this explicitly we write the AB tunneling amplitude for the spin tun- 
neling in the spin coherent state representation: 



icf B Q,r) = <n f \e-^\n l > 

= J V[^T')]V[cos9(r')]e- SE ^ T '^ T ^, (30) 



where 



r/2 

S E W),<f>(?)\= J [£N(0,<f>)-iS<j>(cos9-l)]dT' (31) 

-r/2 



is the Euclidean action including the Berry-phase term. Here, the boundary 
conditions are | >=| >=| ^ ,0 >, | O(0 >=| Cl f >=| 9 ,2tt/N >. 

The analogy of the spin tunneling with the tunneling in S 1 geometry can be 
revealed by making the action (31) similar to Eq. (2), which can be done by 
performing the integral over cos# in Eq. (30). Note, however, that the integral 
is not Gaussian for N > 3. Since we are interested in the general feature of the 
effect of the AB phase on the spin tunneling, we look for an analytical result. 
To do this, we assume that A and h are small enough so that an approximation 

sin^ cos 2 9 (32) 

can hold for N > 3. By using this approximation and performing the 9 integral 
we obtain an effective Euclidean action 



'eff — 



r/2 r/2 

/ dr M^I^ + Vn^) +t J drA N (m (33) 



-r/2 



where 



M N {<P) 



2A(i + M^zi cos 7V0)' 
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S 



1 - 




\NS N - 2 
2 



(34) 



cos 



V N (<t>) 



-AS 2 \S N - 2 cosN<i) + 




XNS N - 2 
2 




COS 



Note that the action is complex. Comparing the real part with Eq. (2), Mtv(</>) 
can be understood as the position dependent mass, and hence Vn(4>) become 
the potential that has an iV-fold symmetry and minima at <fi = 0, 2%/N, ■■, 2(N— 
l)7r/2 on the S 1 because <fi is restricted in the range < <fi < 2ir. The imaginary 
part contributes to the AB tunneling amplitude (30) as a phase factor which 
plays the same role as the phases carried by the instanton and anti-instanton: 
For the tunneling from = to = 2n/N, the instanton picks up a phase 
factor $at = Jq W ^ N dip An {<$)-, while the anti-istanton acquires a phase — $jv- 
When these two paths interfere each other and make a closed loop around 
the hard axis along which the magnetic field is applied, the integral in the 
imaginary part become the AB-phase. From this notion the AB-phase can be 
calculated by integrating the imaginary part along the S 1 ; 



where N = 2,3, 4, 6. As we have mentioned below Eq. (14), it can be seen that 
the phase carried by the instanton Q N is equivalent to <3>/iV. 

Substituting this into Eqs. (24) and (26) we observe several interesting features 
for the level splittings due to the spin tunneling in the iV-fold symmetric 
potential. First, we see that the energy levels split from the ground state 
oscillate as the field varies, and thus intersect each other at certain values of 
h as expected. From the degeneracy conditions in Eq. (27), we note that not 
all levels are simultaneously degenerate in the cases of iV = 3,4, 6 although 
there can be more than one double degeneracy for a given value of <£>. The 
N = 2 case, which corresponds to the 2-fold symmetric potential on the S 1 
geometry, is more interesting. In this case we do not need the approximation 
(32) to calculate the 9 integral in Eq. (30) and are able to find exact condition 
for the degeneracy. Using the condition (27) and Eq. (35) with iV = 2 we find 
that the two ground states are degenerate whenever 



which coincides with the previously found result [4] noticing that A + 2C = k\ 
and AC = k 2 . Thus, the field dependent quenching of the spin tunneling in a 2- 




(35) 




(36) 



10 



fold symmetric spin system can be understood from the viewpoint of Ahronov- 
Bohm effect. 

Second, since we included the Berry phase term in Eq. (31) the AB-phase 
in Eq. (35) provides more information about the degeneracy. To see this we 
consider the field free case, i.e., h = 0. Then, the AB-phase becomes $ = 2irS. 
Using this and Eq. (26) we find following degeneracy conditions depending on 
the spin: 



when N = even; 



Y double degeneracies for half-integer S 

Y — 1 double degeneracies for integer S 

double degeneracies for half-integer and integer S, when N = odd. 



Notice that, for N > 3, levels can be doubly degenerate for both half-integer 
and integer spin and the number of double degeneracies increases as the sym- 
metry increases. When N = 2, the above results show that there exists one 
double degeneracy for half-integer spin. This means that, since there are two 
degenerate ground states in this case, the level splitting due to the spin tun- 
neling disappears for half-integer spin. Thus, the spin parity effect (i.e., the 
Kramers degeneracy) discussed in reference [3] can be explained by the argu- 
ment based on the Aharonov-Bohm effect. 



4 Summary 



We have studied the Aharonov-Bohm effect on the tunneling in an TV-fold 
symmetric potential on the S* 1 geometry. Using the semiclassical methods we 
found that the low-lying levels split from the iV-fold degenerate ground state 
oscillate as a function of the Aharonov-Bohm phase, from which general de- 
genracy conditions depending on the phase are also obtained. By applying 
these observations to the spin tunneling in the spin system which has an N- 
fold rotational symmetry around the hard axis we found oscillating tunnel 
splittings with varying magnetic field applied along the hard axis and degen- 
eracy conditions which are dependent on the spin. From this, the quenching 
of the spin tunneling in the 2-fold symmetric spin system can be interpreted 
as the Aharonov-Bohm effect. 
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